Exact dynamics of finite Glauber-Fock photonic lattices 
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The dynamics of Glauber-Fock lattice of size N is given through exact diagonalization of the 
corresponding Hamiltonian; the spectra {Xk} is given as the roots of the AT-th Hermite polynomial, 
fliv(Afc/\/2) = 0, and the eigenstates are given in terms of Hermite polynomials evaluated at these 
roots. The exact dynamics is used to study coherent phenomena in discrete lattices. Due to the 
symmetry and spacing of the eigenvalues {A^}, oscillatory behavior with highly localized spectra, 
that is, near complete revivals of the photon number and partial recovery of the initial state at given 
waveguides, is predicted. 



I. INTRODUCTION 

Waveguide lattices, that is, arrays of single-mode 
waveguides coupled by evanescent fields, have been the 
focus of considerable interest due to their ability to simu- 
late a variety of quantum effects under negligible decoher- 
ence. Examples of such quantum effects arc Bloch oscil- 
lations in lattices with linearly varying on-site refraction 
index P~H1] , Zeno effect due to a defect in the coupling be- 
tween the first two waveguide in a lattice with otherwise 
homogeneously coupled components [5], random walks 
in homogeneous lattices [SHE], Anderson localization in 
homogeneous lattices where controlled disorder has been 
added [TO] and the effect of isolated defects on quan- 
tum correlations In particular, lattices where the 
coupling is homogeneous are well understood and their 
analytical closed form time evolution is well known [51 1121 - 

E|. 

Recently, a photonic waveguide lattice where the cou- 
pling between adjacent waveguides varies as the square 
root of their position in the lattice has been proposed. 
The propagation of a classical field in a semi-infinite ar- 
ray of this type, the so-called Glauber-Fock photonic lat- 
tice, has been solved in close analytical form by creatively 
mapping the j-th waveguide to the j-th Fock state, this 
has been shown to produce a classical analogue to quan- 
tum coherent and displaced Fock states at the lattice 
output |15j . The quantum correlations of non-classical 
light input have also been analyzed by numerical diag- 
onalization and classical experimental results for single 
waveguide input have been presented for a lattice com- 
posed of sixty waveguides [TO] . 

Here, in Section [IT] it is shown that the finite Glauber- 
Fock Hamiltonian describing an array of identical waveg- 
uides where nearest neighbor couplings varies as the 
square root of the position of the waveguide is exact diag- 
onalizable. The exact result is given in terms of Hermite 
polynomials evaluated at the roots of the iV-th Hermite 
polynomial, where N is the number of waveguides in the 
system. In Section [Til] it is shown that, for a Fock state 
coupled to the zeroth waveguide, there exists an almost 
complete revival of the probability to find the photons 
back in the starting waveguide; the opposite occurs when 
the state couples to the end waveguide, where oscilla- 
tions are polychromatic and revivals are weak. Single- 
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Figure 1. (Color online) Some systems modeled by the so- 
called Glauber-Fock Hamiltonian, Eq.|l]).(a) Array of identi- 
cal photonic waveguides, (b) micro-ring resonators, (c) cou- 
pled cavity arrays, and (d) capacitive coupled strip-line res- 
onator arrays. 



waveguide revivals do not occur in semi-infinite Glauber- 
Fock lattices, nor in uniform lattices unless multi-input 
phenomena or tunning of the lattice is used; say, Talbot 
effect [17] or Bloch oscillations [T][2], in that order. Single 
input revivals occur in more complex waveguide lattices, 
for example, harmonic oscillator [18] . Jaynes-Cummings 
[TO] and Glauber-Fock oscillator [5D] lattices. The dy- 
namics of initial states involving multiple waveguides is 
presented in Section [Tvl Time evolution of two- waveguide 
input, in particular product and NOON states, is ex- 
plicitly discussed and revivals for fidelities, in the single- 
photon superposition case, and two-photon correlations, 
for the two-photon case, are shown. Finally, in Section 
jV] conclusions are presented. 



II. EXACT DYNAMICS 

The Hamiltonian describing a one-dimensional chain of 
N cavities where nearest neighbors are coupled according 
to their position in the chain is given by, in units of h, 



N-l 



JV-2 



H = u + 9 VT+T (a]a j+1 + aja] +1 J (1) 
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The operator at (a^) creates (annihilates) a photon in the 
fc-th cavity, the constants lo and g are the field frequency 
and the base coupling between cavities; these terms are 
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related to the refraction index of the waveguides and the 
inter-waveguide distance in the photonic lattice. Figure[l] 
shows a sampler of systems modeled by this Hamiltonian. 

In the frame defined by the free field, Uo(t) = 
e -iutJ2 j=0 a i a i ^ dynamics is given by the Hamilto- 
nian, in units of kg, 



Hr = 



N-2 

E 

j=0 



VT+T (a]a j+ i + %at +1 J 



(2) 



which Heisenberg equations of motion, dj 
those of a Glauber-Fock photonic lattice 
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In matrix form, the Heisenberg set is given by the expres- 
sion d t a = —iMa where the matrix elements of M are 
m j:k = 5j + i tk Vk + S Jik+1 ^/j, with j, k = 0, 1, . . . , N - 1. 
As the matrix M is real, symmetric and tridiagonal, there 
exists a solution for this differential set given by, 
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where A is a diagonal matrix with elements Xj 
being the eigenvalues of M, and the j-th eigenvector with 
fc-th element Vj k defines the j-th row of the matrix V. As 

v j,k " " " 



Vk.j, it is possible to write 
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aj{t) = V U jik (t)a k {0), 



fc=0 
N-l 

E 

1=0 



e tXlt vijV lk . 



(8) 
(9) 



The collective modes b — Va allow rewriting the equation 
of motion for the annihilation operator, Eq.Q, as 
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which leads to the diagonal Hamiltonian, in units of hg, 

N-l 

Hi =J2 X i%- ( n ) 

3=0 

Thus, the solution to Heisenberg equations of motion lead 
to the diagonalization of the Hamiltonian [21] . 

The method of minors [22] delivers the eigenvalues Xj 
as the roots of the characteristic polynomial of the matrix 
M given by 



H 
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V2 



= 0, j=0,...,A-l. 



(12) 



where H n {x) is the n-th Hermite polynomial [23 which 
zeros are well known |24j . The eigenvectors are given by 



solving the eigenvalue equation (M — IXj)vj — leading 
to the recurrence relations, 

- XjVjfl + Vlv jA = 0, (13) 
vkvj^-i - XjVj ik + Vk + lvj >k +i = 0, (14) 

for j = 0, . . . , N - 1 and k = 1, 2, . . . , N - 2. These 
recurrence relations are fulfilled by the matrix elements 
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Z^k=0 a j,k 



(15) 



(16) 



for j, k = 0, . . . , JV — 1. The Hamiltonian has been diag- 
onalized by giving an analytical closed form. 



III. SINGLE- WAVEGUIDE INPUT 

Let us consider as initial state a Fock state with m > 1 
photons in the p-th waveguide, 



1^(0)) 



1 



i m (0)\Q), p = 0,...,N-l. (17) 



The time evolution of the photon number at the q-th. 
waveguide for such single input initial state is 

(n q ) p ee (v>p(o)|aJ(t)a,(t)|Vp(o)) = ™\u P , q (t)\ 2 - (is) 

For a single-photon coupled to a single- waveguide, it is 
possible to compare published results for a semi-infinite 
lattice [HI [16] with the results presented here for a finite 
large lattice where propagation time is short enough to 
guarantee that the probability of finding the photon near 
the (N — l)-th waveguide is negligible. This leads to the 
following identities, 



it4-,„i 2 = i e - f2/ V)V ( ^4i(Oi 2 <i9) 
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where the expression (x) is a generalized Laguerre 
polynomial [23] . In the case of the initial photon imping- 
ing the zeroth waveguide, the value becomes 
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Figure [2] shows the time evolution of the mean photon 
number for a lattice composed of two hundred waveg- 
uides. The cases of a single-photon starting in the zeroth, 
Fig.([2|(a-c), and fifth, Fig.([2])(d-f), waveguide are pre- 
sented; the figure is equivalent to those found in [T51I16] . 
Here, the focus will be the study of finite size behav- 



ior. Equation ( 18 ) allows the emergence of well defined 
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Figure 2. (Color online) Exact time evolution of the mean 
photon number at the j-th waveguide for an initial state con- 
sisting of a single-photon in the (a-c) fifth and (b-f) zeroth 
waveguide for a lattice size of two hundred waveguides. Snap- 
shots at times (b,e) t = 3 and (c,f) t — 6 comparing the exact 
finite result (blue dots) with the semi-infinite result (solid red 
line). Time is given in units of g . 



oscillations and almost complete revivals at the zeroth 
waveguide as {no)o = "i|£^o.o| 2 with 
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U Q ,o = ^2 w 0: 



ie 



-i\,t 



™o,i 



2=0 



l^k=0 a l,k 
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The weights w 0i ; are highly centralized; that is, revivals 
at the zeroth waveguide are due to a few-component 
chromatic oscillation of the mean photon number; see 
Fig. fcfy (b) . Meanwhile, when the Fock state couples to 
the (N — l)-th waveguide, 
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all eigenvalues contribute equally to the time evolution of 
the mean photon number (njv_i)jv_i = ttj,|J7jv— 1 1 2 
for any given size of the lattice and major revivals have a 
lower probability to appear than in the opposite extreme 
case, (n )o; see Fig. ^ (c) . 

Figure |3^a) shows the time evolution of the mean pho- 
ton number at the j-th waveguide for the N possible ini- 
tial conditions where a single-photon starts at the same 
waveguide, (rij)j = m\Ujj\ 2 . Figure |3|b) and (c) shows 
the spectral weights corresponding to the first and last 
three waveguides, in that order. It is possible to see that 
major revivals of the mean photon number at the initial 
waveguide only occur when the input field impinges near 
the zeroth waveguide. 



Figure 3. (Color online) Summary of the N different cases 
where the initial state is given by a single-photon in the j-th 
waveguide, (a) Only the time evolution of the mean photon 
number at the j-th waveguide, {nj(t))j, is shown. Spectral 
decomposition Wj t k = |(0|aj&|.|0)| for the (b) first j = 1 (black 
circles), 2 (blue squares), 3 (red diamond), and (c) last, j = 
7 (black circles), 8 (blue squares), 9 (red diamond), three 
waveguides. Time is given in units of g~ . 
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Figure 4. (Color online) Exact time evolution of the mean 
photon number at the j-th waveguide for an initial state con- 
sisting of a m-photon Fock state in the zeroth waveguide, noj 
from Eq.p5|, for a lattice size of (a) N = 11 and (b) N = 20 
waveguides. Time is given in units of g . 



Figure [4] shows the time evolution of the mean photon 
number, Eq.([T8]), in a finite Glauber- Fock lattice com- 
posed by eleven, Fi|Qa), and twenty, FijQb), waveg- 
uides for the initial state given in Eq.([T7| at the zeroth 
waveguide, p = 0. In this figure, it is possible to see an 
oscillator-like behavior for single-waveguide input with 
almost complete revivals of the mean photon number at 
the zeroth waveguide; well defined minor revivals occur 
at the last waveguide of the lattice too. 



IV. MULTI- WAVEGUIDE INPUT 

The focus of this section will be the study of 
the oscillator-like behavior, already found for single- 
waveguide input, in the case of multiple- waveguide input. 
The dynamics for single-photon and a pair of relevant 
multi-photon, that is, product and NOON, states will be 
discussed. 



4 




Figure 5. (Color online) Time evolution of the mean photon 
number, Eq.(27l, for a Bell state, Eq. ( 26 1 with a = /3 = 
l/\/2, coupled to the (a)(j,fc) = (0,1) and (b)(j,fc) = (0,4) 
waveguides. Time is given in units of g" 1 . 



A. Single-photon superposition states 

The general state describing the superposition of one 
photon coupled to the lattice is 



Mo), 



N-l 
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1. 



(24) 



Such an initial state covers, for example, Bell states, 
when a single-photon is evenly coupled to two waveguides 
only, and W-states, when the photon is evenly coupled 
to all waveguides. The mean photon number at the g-th 
waveguide evolution for such class of states is given by 
the expression, 



i i)ps = i E ^j.gWI 5 



(25) 



The simplest of such states is the output of a beam split- 
ter with transmission (reflection) a, [fi = yl — a 2 ) cou- 
pled to two waveguides, 



1^(0), 



(aa] + /3al)\0), 



l q ) s = \aU j>q (t) 



(26) 

PU k , g (t)\ 2 . (27) 

Figure[5ja) shows the evolution of a 50/50 beams splitter, 
that is, a Bell, initial state coupled to the (j, k) — (0, 1) 
waveguides and Fig. J5^b) for (j,k) = (0,4). The 
oscillator-like behavior of the mean photon number evo- 
lution is well defined for the first set of initial conditions 
and noisy for the second. The revivals of the mean pho- 
ton number at the starting waveguide are well defined 
and almost complete for both sets as shown in the fol- 
lowing. 

In order to quantify the likeliness between the initial 
state and its time evolution, it is possible to define a 
fidelity, 



j?® = \(ip(o)s\ms 

= \\a\ 2 U^ + 1/ 



2 U k , k + 2Re(a*/3)U. 



3,k\ 



(28) 



that is, the fidelity reaches a value of one if the evolved 
state and the initial state are identical. Building upon 
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Figure 6. (Color online) Time evolution of the fidelity, 
Eq.(28l, coupled to the (j,k) = (0,1) waveguides with (a) 
N = 20 and (b) N = 200 and to the (j, k) = (0, 4) waveg- 
uides with (c) N = 20 and (d) N = 200. Absolute value of 
the normal mode spectra for (e) iV = 20 and (f) iV = 200. 
Time is given in units of <? _1 . 



the spectral decomposition ideas shown in Fig.Q, it is 
possible to realize that, in order to get strong revivals, the 
initial state should be coupled to the zeroth waveguide 
and the coupled waveguides should be as close as possi- 
ble. A lattice of larger size will present stronger revivals 
of the fidelity because the dominant eigenvalues of the 
spectra of Ujj and U k ,k will be closer to each other al- 
lowing them to interfere constructively; of course, period 
of the revivals will increase with the size of the lattice 
and, for real world systems, losses will accumulate. Fig- 
ure [6^a,c) shows the time evolution of the fidelity, for the 
cases mentioned above, Fig{5j and Fig.([6|(b,d) for iden- 
tical initial conditions but a lattice of size two hundred 
waveguides, iV = 200. Figure [6^e,f ) show the absolute 
value of the normal mode spectra of the lattice, {Aj}, 
ordered such that Ao < Ai < . . . < Xjsr-%. 



B. Product states 

Published results for the semi-infinite Glauber-Fock 
lattice include the evolution of the two-photon separa- 
ble state [T5] belonging to the class, 



i<Mo)>=n a ki°>' 



(29) 



where x = (xi, . . . , x k ) with x,; € [0, N — 1] and Xj xj 
for any i =^ j. It is simple to calculate the evolution of 
the photon number at the q-th waveguide, 



k 

E 

3=1 



U x 



(30) 



That is, the probability of finding almost complete re- 
vivals in the oscillation of the mean photon number at 
the starting waveguides depends on the functions U x . q 
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photon number is identical with that of the two-photon 
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Figure 7. (Color online) Time evolution of the two photon 
correlation, Eq. ( 33 1 , for a two-photon product state, Eq. (31 1, 
in a lattice of size N = 20, with (j, k) — (1,2) at times (a) 
t = Ox, (b) 3x, (c) 5x, (d) 6x, (e) 7x, (f) 10 x tt/(10|Aio|)- 
Time in units of 



0.04. 




Figure 8. (Color online) Time evolution of the two photon 
correlation, Eq.( |36[ ), for a two-photon NOON state, Eq.| |34[ ) 
with <f> = 0, in a lattice of size N = 20 with (J, k) — (1, 2) 
at times (a) t = Ox, (b) 3x, (c) 5x, (d) 6x, (e) 7x, (f) 
10 x 7r/(10|Aio|). Time in units of g' 1 . 



having almost identical centralized spectral distributions; 
that is, the criteria commented above regarding dis- 
tance to the zeroth waveguide, separation between field 
mode components, and size of the lattice hold. For two- 
photon product states, the mean photon number evo- 
lution and the two-photon correlation function are well 
known [T5] , 

MOW = a]a\ |0), (31) 
K) ps = \U M \ 2 + |t/ M | 2 , (32) 
r& = \U P jU 9tk + U Plk U q j\ 2 . (33) 

The two-photon correlation can be used to verify the par- 
tial recovery of the initial state. Figure [7] shows the time 
evolution of the two photon correlation for a two-photon 
product state in a lattice of size twenty waveguides. It is 
possible to see the partial recovery at time approximately 
equal to t — 7r/A m i n as expected from all previous argu- 
ments; A m i n = | Ajv/2 I f° r both even and odd (with N/2 
rounded to the next integer) lattice parameter TV and, 
again, Ao < Ai < . . . < Aat-i. 

C. NOON states 

Higher order NOON states are a highly entangled class 
of states, 

|^(0)) = ^= (fit™ + e im *a\ m ) |0) , m = 2, 3, . (34) 

(n q ) = j(\U j , q \ 2 + \Uk, q \ 2 )- (35) 

Notice that the restriction m > 2 has been implemented 
to separate the single-photon superposition treated be- 
fore which delivers an interferometer-like mean photon 
number evolution. It is well known that for the two- 
photon case, m — 2, the time evolution of the mean 



separable state considered above but the two-photon cor- 
relation is different [S] 

r P9 = KjU^l 2 + |t/ p , fe CV| 2 + 2Re (e^U^U^Up^U^f) 

Figure [8] shows the time evolution of the two photon 
correlation for a two-photon NOON state in a lattice 
composed by twenty waveguides. Again, it is possible 
to see the partial recovery at time approximately equal 
to t = 7r/A m i„ as defined above. 



V. CONCLUSIONS 

The exact dynamics of a Hamiltonian describing a fi- 
nite array of coupled identical photonic waveguides where 
the coupling varies as the square root of the position of 
the waveguide in the lattice was presented. 

It was shown that the closed form analytical time evo- 
lution predicts a strong oscillator-like behavior of the lat- 
tice for single and double waveguide input. That is, the 
initial state is partially reconstructed periodically, lead- 
ing to revivals of the mean photon number evolution at 
the waveguides where the photons started. The strength 
of the reconstructions, thus of the revivals, is a combina- 
tion of input distance from the zeroth waveguide, inter- 
waveguide input distance, and lattice size. 
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